Stability under Perturbations of the Large Time Average Motion of Dynamical 
Systems with Conserved Phase Space Volume. 
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The stability against perturbations of a dynamical system conserving a generalized phase-space 
volume is studied by exploiting the similarity between statistical physics formalism and that of 
ergodic theory. A general continuity theorem is proved. Resulting from this theorem the double 
average - time and ensemble - of an observable of a weakly perturbed ergodic dynamical system 
is only slightly changed, even in the infinite time limit. Consequences of this statistical analogue 
of the structural stability are: extension of the range of practical applicability of the Boltzmann 
Ergodic Hypothesis, justification of the perturbation method in statistical physics, justification of 
the numerical approximations in molecular dynamic calculations and smoothness of the transition 
from bounded to unbounded motion as observed in numerical simulation of anomalous transport in 
tokamaks. 
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The stability of the predictions from equilibrium sta- 
tistical physics is closely related on one hand to the Botz- 
mann Ergodic Hypothesis (BEH) 0] and on the other 
hand to the persistence of the ergodic properties when 
a dynamical system (DS) is exposed to small perturba- 
tions. The accompanying continuity problem is related 
to the accuracy of the BEH in the neighborhood of a DS 
with proven ergodicity Q 0. 

There are known examples of DS's, for which the ra- 
tional/irrational character of the parameters leads to a 
change in the behaviour from non chaotic to mixing j^. 
Weakly chaotic DS 0, or DS at the verge of chaoticity 
0, are important for the study of anomalous transport 
in tokamaks. A basic problem related to perturbations 
is to guarantee the validity of the numerical approxima- 
tions used to study the DS or, more generally, to prove 
the continuity of the infinite time limit with experimen- 
tally uncertain parameters. Such problem arises for ex- 
ample in area preserving maps {e.g. for the standard 
map 0) where the area of the phase space visited by the 
trajectories starting from a fixed small domain is a dis- 
continuous function of the parameters. This behavior is 
the consequence of the destruction of invariant circles at 
well defined critical values of the stochasticity parameter. 

We shall prove, using a main theorem, that a slightly 
perturbed DS defined in a finite phase space volume, de- 
spite not necessarily being ergodic, behaves almost like 
an ergodic DS when probed with observable relevant from 
the point of view of statistical physics. For the DS de- 
fined in a infinite phase space volume as in the case of 
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the dissipative DS _£] (i. e. when the trajectories are 
unbounded with 100% probability as for example in the 
case of the infinite ergodic DS's), the theorem justifies 
the absence of threshold in the transition from bounded 
to unbounded motion, as observed in numerical study of 
anomalous particle transport in tokamak . 

The continuity of the infinite time limit of the double, 
time and ensemble averages (ITEA) is studied. It will 
be proved that this double average is stable under small 
perturbations in the ergodic and the dissipative cases. 
For DS that are not ergodic or dissipative, it is possible to 
have discontinuous behaviour, as in the previous example 
of the standard map. As a consequence of our analysis, 
unphysical number theoretic aspects of ergodic properties 
of weakly stochastic systems are clarified. The continuity 
results proved below is an additional argument in support 
to perturbation method in equilibrium statistical physics 

and prove the stability of the results of microcanonical 
molecular dynamics simulations when the unperturbed 
system is ergodic. 

For the sake of simplicity and generality, we consider 
an abstract, discrete time, measure preserving, dynami- 
cal system Ij . Such system may be defined on the phase 
space M with the generalized phase space volume fi as 
conserved measure, and a discrete time one to one evolu- 
tion given by a (stroboscopic) map x — > r(x). Typically, 
X = (Piq)7 and d^(x) = dpdq. For the sake of mathe- 
matical rigor, we denote by A the family of the subsets 
of M, where /i is defined. For every subset A of the fam- 
ily A, we postulate that ^jl{A) = ^[T{A)] 0, which is 
nothing else than the well known Liouville theorem in 
the case of Hamiltonian systems. We shall denote the 
DS by {AI, A, fJ,,T). The measure can be either finite 
when it is normalized to /i(Af) — 1, as in the case of the 
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classical finite spatial extent DS with a potential energy 
bounded from below or infinite like in the case of a 
two component classical plasma 0] . The definition of an 
ergodic DS in the mathematical literature, for both the 
finite and the infinite cases, is the following: if a subset 
A of M is invariant, then either A or its complement has 
zero measure. From this definition follows the equality 
of time and ensemble averages Q, 0- 

For any subset A of A, we denote by v{T, A) the subset 
of M visited by the trajectories of the map T, starting 
from A. The continuity of /x[z/(T, A)] as a consequence of 
the main theorem will now be considered. 

Without loss of generality, the relevant properties 
of the DS's from the statistical physics point of view 
can be studied, in a technically convenient manner, in 
Hilbert space H of square integrable functions L'^{M,^) 
with the scalar product defined as (</> | V') = 
^*(x) ^(x) d^(x). The hnear evolution operator U 
(also denoted U{T)) acting in H is defined by ?/'(x) 
-0 [T(x)] = {U '4'){'^)- It is unitary due to the conser- 
vation of the measure /z. Choosing (/'(x) as the ini- 
tial time probability density function (PDF) po(x), and 
-0(x) = r(x) an observable, then (0 1 1/* V) = {po I Y) 
is the ensemble average of the observable Y{yi) at time t. 
Consequently, the physical meaning of the mathemat- 
ical concept of weak convergence of DS can be made 
obvious. According to Refs. 0, the parametrized 
family of DS's (Af, ^, ^,Te) converges weakly to the DS 
{M,A,y,,T) iff for any 0(x),V'(x) belonging to TL it fol- 
lows that {(j) I U{T^) I U{T) i)), when e ^ 0. This 
weak convergence is denoted -> T. We note that the 
weak convergence follows from the usual convergence of 
the maps x — > Te(x), but the converse, in general, is not 
true. 

The elements of the one parameter family of DS's 
(Af, /i,Te) can be considered as weak perturbations 
of the DS {M, A, fi, T). This is the case, for instance, for 
hard ball systems obtained by adding a small potential 
energy or by inserting small reflecting obstacles of zero 
volume. It is easy to prove, using simple Hilbert space 
algebra, that Ti*"' ^ T^^^ is a consequence of ^ T, 
where T'^^'^ is the k — th iterate of the map T. Accord- 
ing to the definition and the physical interpretation of 
the weak convergence, it follows that for any abstract 
dynamical systems the time evolution of the ensemble 
average at any finite time has a continuous dependence 
on the weak perturbation of the unit time evolution. This 
result being established, it remains to study the effect of 
a weak perturbation on the statistical properties of the 
infinite time limit. 

The infinite time limit (iV — > oo) of the ensemble aver- 
age (po \U^Y)^ J^^ po(x) Y [r(^)(x)] d^i{^) in general 
does not exist at the exception of mixing dynamical sys- 
tems jH- Consequently we must restrict ourselves to the 
weaker information given by the ITEA, which is typically 
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Recalling the von Neumann ergodic theorem this limit 
exists if 0(x) and ipi^) are square integrable 0. 

Moreover, denoting by Hinv the subspace of Ti. con- 
sisting of invariant functions, we have that tpi^) belongs 
to H,nv iff ^/'(x) = [T(x)] (i.e. iff V = U{T)i}). In a 
similar way, denoting by the orthogonal complement 
of Tiinv j-e. Ti ^ Ti-inv © Ti-c and defining the operator 
P{T) as the orthogonal projector on the subspace TLmv, 
the von Neumann ergodic theorem states: 
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In particular, if the initial time PDF is po(x) then the 
function /3oo(x) = [P(T)po](x) is the PDF of the points 
when the probability of a domain is calculated according 
to the mean visiting time fraction. The problem of the 
continuous dependence of the ITEA, in the sense of the 
weak convergence, reduces to the study of the following 
question: what are the conditions under which Tg ^ T 
when £ ^ leads to {(j) \ P{T,) ip) ~* P{T) ip) for all 
</>(x),i/;(x) belonging to H, i.e. P{T,) ^ P{T). 

A theorem is the starting point for a unified study of 
the finite or infinite DS. 

Theorem 1: Consider the sequence of DS's 
(Af, A, /i, Tg) that are weak perturbations of (Af, A, ji, T), 
i.e. suppose that —> T. Then for all </?(x) be- 
longing to Ti. and V'c(x) belonging to Tic, we have 

(^|[p(Te)-p(r)]Vc) ->o. 

The proof is performed in three steps: 

A) The statement T is equivalent to 
{ipi I [C/e - U] ip2) — > for any (pi,(p2 belonging to H. 
(We use here the shorter notations Ue — U{T^), U{T) = 
U, P{T^) = P, and P(T) = P). Since Ue and U are uni- 
tary, this statement is also equivalent to an apparently 
stronger statement: for any ip{'sc) belonging to TL we have 
\\{Ue ~ U)ii}\\ ^ for e ^ 0. 

B) We denote 77(e) = \{ip\{Pe ~ P)^c)\, then for e 
Tic i.e. Pipe = we have 77(e) — \{lp\Ps ipc)]- We now 
consider the normalization \\ip\\ = 1. The key lemma used 
to prove the von Neumann ergodic theorem 'I:] states that 
for ipc belonging to Tic there exist a xs ^ such that 
llV'c -{U -I) X5 II < 5/2. We obtain 

r7(e) < \{^\P,[^,^{U -1)X8])\ 

+ \{PM[{U-1)-{U,-1)]X8)\ 
+ \{^\P,{U,-\)X5)\ ■ 

Using the relation P^ = Pe iU the last term is zero. 
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C) From Schwartz inequality it follows that 77(e) < 
(5/2+ ||(J7e — C/)X(5|| ■ Indeed, from A) there exist a e such 
that \\{U^ - U)xs\\ < S/2. Then, for 6 arbitrary small 
there exist a e such that \{^p\{P^ — P)^pc)\ — r]{e) < 6. 
This completes the proof. 

We are now in position to prove a second Theorem 
valid for abstract ergodic systems. 

Theorem 2: Consider the DS's {M,A, IJ.,T^), 
{M,A,n,T), ^ T, and suppose that the DS 

{M,A,iJ,,T) is either ergodic (when ^{M) — 1), or 
dissipative (when ^{M) — 00). Then for all (^(x) and 
'0(x) belonging to Ti we have: 

(V? I [P(Te) - F(r)]?A) ^ fore 0. (2) 

Due to the linearity of the scalar product and the pre- 
viously used Lemma it is sufficient to get a proof for 
the case when 'ip(x.) = ^/;j„„(x) is an invariant func- 
tion. In our case ipinvi'x.) = const Q. We then have 
U{Tf)'4>inv = U{T)ipinv = Anv and by using (P) it fol- 
lows that P{T^) ipinv — P{T)'ipinv — 'ipinv, which com- 
pletes the proof. 

According to the standard terminology of statistical 
physics, a DS associated to an autonomous Hamilto- 
nian systems with N degree of freedom, is considered 
ergodic if the restriction of the motion to the hypersur- 
faces determined by a fixed energy is ergodic in mathe- 
matical terms. We denote by -ff(x) and H^(x.) the un- 
perturbed (assumed ergodic) and the perturbed Hamil- 
tonian functions, respectively. Then ipinvi'x.) is of the 
form ?/'int)(x) = /(i?(x)). Without loss of generality, we 
consider the smooth functions ipinv{^) and denote the 
quantity Je(x) = /(iJ(x))//(i?e(x)). Then, if in addi- 
tion to the conditions of the Theorem 1 we assume that 
Je(x) — 1| < S{e)/4: with lime_^o'5(£) — (this happens, 
for instance, when f{H{x.)) = exp(— /3 H{x.))), we obtain 
the physical result that: 

An ergodic (in physical terms) Hamiltonian system sat- 
isfying the above conditions on Je(x) obeys J^). 

Indeed, recalling Theorem 1, we only need to compute 
P{T^) 4'inv and by the definition of P{T^) = P^ given by 
n we must estimate (C/^V-'mi>)(x). Simple algebra then 
leads to iU^4>^nv)i^) = V^,™(x) Jari'-')(x))/jax) = 
■0i™(x)Le,fc(x). We also have |ie,fc(x) - 1| < 6(e), for 
small (5(e). From the normalization ||iy9|| = HV'inull — 1 
it follows that |((y5 | [U^ — U'^] ipinv)\ < <>{e) and conse- 
quently that \{ip I [P(Te) - P{T)] 'ijj,„v)\ < '5(e). This re- 
sult completes the proof. 

We shall denote in the following by 1a (x) the func- 
tion equal to 1 in the domain A but zero outside (the 
characteristic function of A). If the system is ergodic, 
then P(T) 1a(x) = ^(A) 1a./(x) = n{A). As an applica- 
tion, we consider fi{M) = 1 and (p{x) = ls{x)/^{S) the 
normalized distribution function of the starting points in 



the previous Eq. For simplicity, we consider the 

"source " of the particles localized in the domain S. We 
assume there is a "detector" in the domain D, measuring 
the mean visit time fraction, spent in D. According to^ 
this quantity is given, in general, by {ip | P(re) l£)(x)), 
which is close to {^p \ P{T) l£)(x)), or, due to the ergod- 
icity of T, close to ^Ji{D). Consequently, irrespectively of 
the positions of the "source" and of the "detector" and 
for sufficiently small perturbations the counting is almost 
the same as that performed on the unperturbed ergodic 
case. There are no possibility to protect the detector 
from the particles coming from the source. // the system 
is not ergodic, then in general this property is no longer 
true, as in the example of the near critical standard map. 
A striking example of non ergodic system having a sen- 
sible dependence on the perturbations is given in 10]. It 
is shown there how analytically calculated perturbations 
can lead to a drastic reduction of the particle transport 
in fusion plasma. 

An alternative way to define ergodicity for finite mea- 
sure DS {p{M) = 1) is to impose fJ,{v{T, A)) = 1 to every 
subset A with ^{A) > 0. This means that the trajectories 
starting from an arbitrary small neighborhood of a phase 
space position completely fill the whole phase space. We 
prove the expected property: the trajectories of the per- 
turbed DS will almost fill the entire phase space. To get 
the proof, we first observe that from Theorem 1 we get a 
result on strong continuity: 

Let {M,A,n,Te), {M,A,fi,T), T, ^ T, and 
{M,A,^,,T) he ergodic, then for any ipi^) belonging to 
H we have \\[P{T,) - P(T)] V'H ^ 0. 

This result is a direct consequence of (O and of the 
fact that P{Te) and P(T) are projectors. 

From this remark, we get the following result on ap- 
proximate ergodicity: 

Let {M,A,tx,T,), {M,A,ii,T), T, ^ T, /x(M) = 1, 
and (M, A, T) ergodic. Consider a subset S such that 
fi{S) > 0. Then fi{v{T^, S)) -> 1. 

This proof follows from the use the strong continuity, 
with = ls(x), and from the fact that P(T£)l5 is 

zero outside the domain v{Tg,S) and that P(T) l5(x) 
is a nonzero constant function. This result shows that 
for a slightly perturbed ergodic DS the trajectories fill 
almost the same phase space as the trajectories would do 
in the case of the unperturbed DS. Again, if the system 
is not ergodic then, in general, this property is no longer 
true. A counterexample is again the standard map at the 
critical value of the stochasticity parameter. 

For clarity, let the phase space M be a torus, with 
coordinates {p, q) defined modulo 1 and define a map 
{P, q) T{p, q) by p ^ p' = p + a, q ^ q' = q + P 
(mod 1). When a, 0, and a/P are all irrational, then 
the map is ergodic p]. Let the perturbed sequence T^ 
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given by the non ergodic, rational approximations of a w 
Qfj = N^/D^, with irreducible fractions. If we choose 
two small circles of diameter i5 for the source S and the 
detector D, then if > 1/5, we get v{T^,A) = M, so 
that the trajectories starting from S visit every parts of 
the phase space. From this point of view, the system 
behaves like an ergodic one, but the frequency of visit 
of the detector approaches a limiting value only when 
^ 1/6. This example illustrates the fact that despite 
a weak perturbation of an ergodic DS does not necessarily 
lead to an ergodic DS, the latter behaves like an ergodic 
one. 

The previous study of the very simple dynamics on the 
torus is of special interest because the dynamics associ- 
ated to quasiperiodic time dependent Hamiltonian, with 
f degrees of freedom, having N frequencies, can be re- 
formulated as an autonomous DS in an enlarged phase 
space, with 2f + N dimensions, and N cyclic coordi- 
nates: X — {pi, ...,pf,qi, ...,qf, (pi, (^jv} and the in- 
variant measure dfj,{x)= dp dq dip. This imbedding, 
known as skew product Jj/, allows for the use of the meth- 
ods of ergodic theory. 

The deduction from Theorem 1 is correct whenever 
the invariant functions of T are constants and can be 
used also in the framework of non equilibrium statisti- 
cal physics. In the case of DS with fJ.{M) = oo, when 
(M, A, /i, T) is dissipative, the single invariant function 
is zero i. e. P{T)^ = 0. Then from T^ ^ T and from 
13 it follows that {4>\ P{Te) ip) — > 0. As an example we 
consider the transversal two dimensional guiding center 
motion of a charged particle in a constant magnetic field 
and in a fluctuating electric field reduced to the model 
proposed in 0. The motion is described by a time peri- 
odic Hamiltonian system with a single degree of freedom. 
The typical phenomenon mentioned there is the transi- 
tion form bounded motion to unbounded motion for large 
wave amplitudes with increasing wave amplitudes. 

Using numerical simulations, "the lack of stochasticity 
threshold and a smooth increase of the measure of the 
chaotic regions" (with increasing wave amplitudes) was 
discovered. To explain this effect, we assume that the tra- 
jectories evolving under the action of T or of T^ start from 
the bounded domain A. We also assume that T is dissipa- 
tive (unbounded trajectories) and T^ — > T . In this case 
the invariant functions are zero and Theorem 2 holds. 
We have P(T)1a = 0, and we assume ad absurdum, that 
the trajectories, under the action of are all contained 
in the same bounded domain B. Then {lM\P{Te)lA) — 
(PiT,)lM\lA) = { 1m\1a) = m(^) ■ But PiT,)lA = 
outside of B, so {lM\PiT^)lA) = {lB\PiTe)lA)- K fol- 
lows that {1b\P{T^)1a) = m(^) > 0. Recalling the 
continuity theorem we get (1b| [P(Te) - P(r)] 1a) = 
(Is I P{Te) 1a) — > 0. From this contradiction it follows 
that the transition to unbounded motion is not caused by 
a sudden opening of some transport barriers that protect 



the domain A. Since the quasiperiodic perturbations can 
be treated within our formalism, it follows that the same 
smooth transition will be observed in the case of electric 
perturbations with several incommensurate frequencies . 

As a consequence of theorem 1, it is proved, for er- 
godic and dissipative dynamic systems exposed to small 
perturbation, that the properties that are relevant in the 
framework of the equilibrium and non equilibrium sta- 
tistical physics are only weakly perturbed. In particular 
we solved the problem of the dependence of the ergodic 
properties on number theoretic properties of the parame- 
ters. This continuity explains some aspects of the empir- 
ically observed behaviour of the ergodic and dissipative 
DS when exposed to small perturbations. If explored 
with coarse grained instruments, the behavior of the per- 
turbed ergodic system is almost ergodic. Together with 
the results from 10], we conjecture that between the two 
extreme cases: from the ergodic to the complete inte- 
grable DS, there is an increased possibility to reduce the 
chaoticity. All of the results are consequences of the 
weakly continuous dependence on the perturbations of 
the projector on the invariant states. 

Discussions with R. Balescu, D. Carati (ULB, Brus- 
sels), J. H. Misguich and J. D. Reuss (DRFC, CEA- 
Cadarache), M. Vittot (CNRS, Luminy) are greatly ac- 
knowledged. 
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